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STAR FORMATION RATE INDICATORS IN THE SLOAN DIGITAL SKY SURVEY 
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ABSTRACT 

The Sloan Digital Sky Survey (SDSS) first data release provides a database of « 106000 unique 
galaxies in the main galaxy sample with measured spectra. A sample of star-forming (SF) galaxies are 
identified from among the 3079 of these having 1.4 GHz luminosities from FIRST, by using optical spectral 
diagnostics. Using 1.4 GHz luminosities as a reference star formation rate (SFR) estimator insensitive 
to obscuration effects, the SFRs derived from the measured SDSS Ha, [On] and u-band luminosities, 
as well as far- infrared luminosities from IRAS, are compared. It is established that straightforward 
corrections for obscuration and aperture effects reliably bring the SDSS emission line and photometric 
SFR estimates into agreement with those at 1.4 GHz, although considerable scatter (« 60%) remains 
in the relations. It thus appears feasible to perform detailed investigations of star formation for large 
and varied samples of SF galaxies through the available spectroscopic and photometric measurements 
from the SDSS. We provide herein exact prescriptions for determining the SFR for SDSS galaxies. The 
expected strong correlation between [On] and Ha line fluxes for SF galaxies is seen, but with a median 
line flux ratio Fipm/Fna = 0.23, about a factor of two smaller than that found in the sample of Kennicutt 
(1992). This correlation, used in deriving the [On] SFRs, is consistent with the luminosity-dependent 
relation found by Jansen et al. (2001). The median obscuration for the SDSS SF systems is found to be 
^4hq = 1.2 mag, while for the radio detected sample the median obscuration is notably higher, 1.6 mag, 
and with a broader distribution. 

Subject headings: catalogs — galaxies: evolution — galaxies: starburst — radio continuum: galaxies 



1. INTRODUCTION 

The current star formation rate (SFR) of a galaxy is 
one of many important parameters used in developing our 
understanding of galaxy evolution. Various different indi- 
cators of galaxy SFR exist at different wavelengths, and 
include Ha and [On] emission line luminosities, ultraviolet 
continuum luminosity, far-infrared (FIR) luminosity and 
radio luminosities (see reviews by Kennicutt 1998; Con- 
don 1992). There are also strong suggestions that X-ray 
luminosity is an important SFR indicator (Grifiths et al. 
1990; White & Ghosh 1998; Ghosh & White 2001; Ptak et 
al. 2001; Brandt et al. 2001; Georgakakis et al. 2003), al- 
though very faint X-ray observations are typically required 
to detect X-ray emission from star formation processes. In 
recent years several comparisons between different star for- 
mation indicators at varying wavelengths have been made, 
primarily establishing broad agreement between each but 
with detailed discrepancies and large scatter in the rela- 
tions (Bell 2003; Buat et al. 2002; Hopkins et al. 2001; 
Sullivan et al. 2001, 2000; Cram et al. 1998). One of the 
limiting factors in these comparisons to date is the mostly 
heterogeneous nature of the data being compared, and the 



relatively small numbers of objects investigated, no more 
than a few hundred. 

The Sloan Digital Sky Survey (SDSS, Fukugita et al 
1996; Gunn et al. 1998; York et al. 2000; Hogg et al. 2001; 
Stoughton et al. 2002; Smith et al. 2002; Pier 2003) elimi- 
nates these limitations. The SDSS is supplying the astro- 
nomical community with images and spectra providing an 
immense resource for use in numerous studies of galaxy, 
quasar, stellar and solar system properties. By identify- 
ing star-forming (SF) galaxies catalogued by the SDSS, a 
very large, homogeneous sample can be used to investi- 
gate the properties of star formation in galaxies. To sup- 
port such studies, we investigate herein the consistency 
of SFRs derived from SDSS Ha and [On] line measure- 
ments with those derived from the obscuration indepen- 
dent 1.4 GHz and FIR luminosities. Further, we empiri- 
cally derive a non-linear calibration of SFR from the SDSS 
u-band luminosity, which gives SFR estimates highly con- 
sistent with the other four estimators. X-ray luminosities 
are not pursued here as an SFR estimator since not enough 
deep X-ray data is presently available for an exploration 
of a homogeneous X-ray detected SF galaxy sample. Re- 
sults from an analysis of ROSAT All Sky Survey (RASS) 
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identifications with SDSS galaxies in the early data re- 
lease (Schulte-Ladbeck et al. 2003) confirm that almost all 
the RASS X-ray sources are identified with known galaxy 
clusters, quasars or active galactic nuclei (AGN). 

In § 2 we describe the details of the current sample. The 
five SFR indicators are each explored in § 3, beginning with 
the 1.4 GHz estimator in §3.1. The Ha SFRs, including 
the obscuration and aperture corrections required, are pre- 
sented in § 3.2. SFRs from [On] , FIR and u-band luminosi- 
ties are each given in §§ 3.3, 3.4, and 3.5 respectively. A 
discussion of the absolute SFR calibrations and the prop- 
erties of the radio detected SF galaxies are explored in § 4. 
We summarise our results in § 5. Throughout this paper 
we assume a (Am = 0.3, 17a = 0.7, H = 70) cosmology. 

2. SAMPLE SELECTION 

The SDSS sample of spectroscopically observed main 
galaxies (Strauss et al. 2002) taken from the first data re- 
lease (DR1, Abazajian et al. 2003) was used as the start- 
ing point for constructing our sample. This catalogue 
of w 106000 unique galaxies yields 3079 galaxies with 
1.4 GHz measurements from the Faint Images of the Ra- 
dio Sky at Twenty centimeters catalogue (FIRST, White 
et al. 1997). The radio identification from FIRST is based 
on positional matching within a radius of l'/5 of the SDSS 
position (Ivezic et al. 2002). These galaxies form our pri- 
mary sample for the investigations presented here. 

Within our primary sample, we made use of three spec- 
tral diagnostic diagrams to classify galaxies as either SF, 
active galactic nuclei (AGN) or unclassified. Following 
Kewley et al. (2001) and Miller et al. (2003), the diagrams 
used were [Om]/H/3 vs [Nil] /Ha, [Om]/H/3 vs [Sn]/Ha, 
and [Oiii]/H/3 vs [Oi]/Ha. We used the conservative re- 
quirements that for a galaxy to be classified as SF or 
AGN it had to be so classified in all three of the diagnos- 
tic diagrams, and remaining unclassified otherwise. Only 
emission lines where the measured flux was greater than 
twice the flux uncertainty (similar to a 2a threshold) were 
considered in producing these classifications (Miller et al. 
2003). Additionally, some objects could also be classified 
as AGN on the basis of having a large ratio of two emis- 
sion lines for a given diagram, even though the other pair 
of lines for that diagram were not measurable (see discus- 
sion in Miller et al. 2003). Saturated emission lines will 
not affect our sample since saturated line parameters are 
not measured, so systems with saturated Ha, for exam- 
ple, can not be classified as SF galaxies. These systems 
are also extremely rare, and mostly dominated by nuclear 
emission from very nearby galaxies, so very few SF galaxies 
are likely to be omitted by excluding such objects. Our 
primary sample consists of 791 SF galaxies, 379 AGNs, 
672 "quiescent" galaxies (with no measurable diagnostic 
emission lines), with the remaining 1237 galaxies either 
inconsistently classified, or lacking the required emission 
lines to be able to be classified in all three diagrams (see 
Figure 1). The ratio of AGN to SF systems is relatively 
high as a result of requiring radio detection (and hence 
a bias towards AGN systems) for inclusion in the sample. 
While a number of the "quiescent" galaxies are likely to be 
absorption line systems, typical of old stellar populations 
in elliptical galaxies, many are likely to be higher redshift 
objects whose emission lines fall below the sensitivity lim- 



its of the SDSS. 

Although true absorption line systems, particular those 
at higher redshifts (z > 0.3), will incorporate the luminous 
red galaxies (LRGs, Eisenstein et al. 2001), our sample has 
been defined to exclude these objects, being drawn from 
the main galaxy sample only (Strauss et al. 2002). This 
point is relevant since LRGs were noted by Ivezic et al. 
(2002) to have a relatively high fraction of radio coun- 
terparts. However, even if the LRGs are not excluded a 
priori, the results presented here do not change at all, since 
almost none of the LRGs satisfy our spectroscopic classi- 
fication to be considered as SF galaxies. 

Of the systems with the necessary emission lines, rela- 
tively few show conflicting classifications, as seen in Fig- 
ure 1. The 177 galaxies inconsistently classified between 
AGN and SF may be composite systems hosting both 
types of activity (Hill et al. 2001), and arc thus inter- 
esting in their own right, although they are not further 
investigated here. 

All the classified SF and AGN galaxies have measurable 
Ha luminosity (£h q ) and Balmer decrements, a result of 
requiring Ha and H/3 for classification, and we restrict the 
current investigation to the 791 galaxies classified as SF in 
all three spectral diagnostic diagrams. To calculate w-band 
luminosities we chose to use the Petrosian tt-band abso- 
lute magnitudes, M u , (calculated using the k-corrections 
of (Blanton et al. 2003)), treating them as AB magnitudes 
(Fukugita et al 1996). Of the 791 SF systems, 752 have 
measurable [On] emission. 

The completeness of the final spectroscopically classi- 
fied sample is complex to define. The initial SDSS spec- 
troscopic sample is complete to an optical flux limit corre- 
sponding to about r = 17.7, and by requiring the presence 
of the [Sli] emission lines for classification, we implicitly 
impose an upper redshift limit of about z — 0.36. By 
also requiring a 1.4 GHz detection from FIRST, we im- 
pose a second flux limit, this time corresponding to about 
0.5 mJy The full distributions in redshift, 1.4 GHz flux 
density and luminosity of the sample are presented in Fig- 
ures 2 to 4. All three Figures show the distributions for 
the whole sample of 3079 galaxies, distinguishing the 379 
AGNs, 791 SF, and 672 "quiescent" galaxies. The pro- 
portion of galaxies classified as "quiescent" becomes more 
dominant as redshift increases. This comes about as a 
combination of the detection limit of the survey and the 
necessary diagnostic emission lines moving out of the ob- 
servable wavelength window. As a result, some galaxies 
classified as "quiescent," especially those at higher red- 
shifts, may not actually be quiescent at all. This is em- 
phasised in Figure 4 where a large fraction of the brightest 
radio sources are classified as "quiescent." In reality these 
systems may possess strong emission features that are be- 
low the detection threshold or lie outside the accessible 
SDSS window. The distribution in 1.4 GHz luminosities 
for the 3079 radio detections spans the very broad range 
from quite faint SF systems, log(Li.4GHz) ~ 20.5, to pow- 
erful radio galaxies, log(Li. 4 GHz) ~ 26 (compare with the 
range found by Sadler et al. 2002, for example). The lu- 
minosity range spanned by the spectroscopically classified 
AGNs is not too dissimilar from that of the SF systems, 
and this is a result of two effects: (1) requiring the pres- 
ence of the necessary emission lines for the spectroscopic 
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diagnostics limits the classified AGNs to lower redshifts, 
and hence lower luminosities; and (2) the more powerful 
AGNs may be dominated by elliptical type host galaxies, 
with absorption spectra. 

Since the results of our SFR analysis do not rely on hav- 
ing a complete sample, and in order to retain the maximum 
number of objects for the analysis, we have chosen not to 
exclude any galaxies other than those required by our se- 
lection criteria above. Although our full spectroscopically 
classified sample thus remains incomplete at the higher 
rcdshift and lower flux density extremes, we can still de- 
fine regimes over which our sample does approach com- 
pleteness. The initial spectroscopic sample is known to be 
highly complete over the redshift range 0.05 < z < 0.1 
(Gomez et al. 2003). The flux density limit to which the 
FIRST survey is complete is about 2mJy (White et al. 
1997), so out to z = 0.1, the 1.4 GHz detections should 
be complete above about log(Li.4GHz/WHz _1 ) = 22.5 
(corresponding, using the calibration described below, to 
SFRi.4GHz ~ 2OM yr~ 1 ). Despite this fairly high lower 
limit for completeness in SFR, the FIRST detections ex- 
tend below the completeness limit to flux densities of about 
0.75 mjy, allowing SFRs of order unity to be probed, al- 
though the radio detected sample at these levels is highly 
incomplete. 

3. STAR FORMATION RATES 

Estimating a galaxy's current SFR is typically done by 
applying a scaling factor to a star formation sensitive lumi- 
nosity measurement for the galaxy (Kennicutt 1998; Con- 
don 1992). Here we use this technique to estimate SFRs 
using Ha, [On], it-band, 1.4 GHz and FIR luminosities. In 
all the SFR calibrations given below we assume a Salpeter 
initial mass function (IMF) and a mass range from 0.1 to 
100 Mq. Changing the mass range or choosing a different 
IMF will, of course, alter the values for the derived SFRs, 
and a discussion of these effects can be found in Kenni- 
cutt (1998). The radio luminosities have been calculated 
assuming a power-law spectral index S v oc appropri- 
ate for galaxies dominated by synchrotron emission, and 
the optical luminosities are calculated by incorporating the 
k-corrections as measured by Blanton et al. (2003). Below 
we describe the details of calculating the SFRs for the five 
estimators we explore, and comparisons between each are 
performed. 

3.1. 1.4 GHz Star formation rates 

Long-wavelength SFR estimates are insensitive to dust 
obscuration, increasing their attraction for SFR investi- 
gations, but they are not without limitations. Radio lu- 
minosity can be generated by AGN as well as star for- 
mation processes, and indeed the majority of apparently 
bright radio sources are AGN. By selecting for star for- 
mation directly from optical spectroscopic features, how- 
ever, we have eliminated this potential source of confusion. 
The detailed physics involved in the connection between 
SF and radio emission, in addition, is still poorly under- 
stood, despite attempts by numerous models to explain it 
(e.g., Condon 1992; Price & Duric 1992; Chi & Wolfendale 
1990). Despite this, radio luminosities, in part because of 
the tight correlation with FIR luminosities (e.g., de Jong 
et al. 1985; Condon et al. 1991; Yun et al. 2001), appear 



to be robust and efficient SFR estimators (Condon 1992; 
Cram et al. 1998; Bell 2003). This, combined with the fact 
that the 1.4 GHz luminosity is insensitive to dust obscu- 
ration makes the use of SFR1.4GHZ very attractive, and 
we adopt this as a reference SFR when investigating the 
details of the other SFR measures. After the detailed in- 
vestigation of the radio-FIR correlation by Bell (2003), we 
adopt the calibration derived therein between 1.4 GHz lu- 
minosity and SFR, 



SFRi.4GHz (Moyr- 1 ) 
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and L c = 6.4 x 10 21 WHz 1 . This calibration produces 
SFRs about a factor of two lower than the calibration of 
Condon (1992) for luminosities above L c , while for fainter 
luminosities the SFRs progressively converge, with the 
present calibration eventually producing larger SFRs than 
that of Condon (1992) below about 3.4 x 10 20 WHz" 1 (see 
discussion in Bell 2003). 

Most of the current generation of sensitive radio surveys, 
including FIRST, are conducted with interferometric tele- 
scopes, which, lacking short-spacing information, have a 
limited sensitivity to emission from more extended struc- 
tures. The extent of this effect is dependent on the short- 
est baseline of the instrument. In particular, the FIRST 
catalogue starts losing sensitivity to 1.4 GHz emission for 
galaxies larger than about 10". For galaxies with a radius 
of 12" FIRST is only sensitive to about 84% of the emis- 
sion, decreasing further for larger sources (Becker et al. 
1995). If this effect is not accounted for, it can result in 
an apparent overestimate of the comparison SFR at low 
SFRs (see, for example, Figure 1 of Cram et al. 1998), 
although it is actually an underestimate in the 1.4 GHz 
derived SFRs. As shown in Figure 5, where the galaxy 
size is shown as a function of SFRi.4GHz from FIRST, this 
starts becoming an issue in the current sample for galaxies 
with SFR1.4GHZ ^ 10M Q yr _1 . To avoid such underesti- 
mates in SFRi 4 GHz for our analysis, we make use of the 
NRAO VLA Sky Survey (NVSS, Condon et al. 1998), a 
1.4 GHz survey made with the VLA in the more compact 
D configuration. As a result of the compact VLA configu- 
ration used for the survey, the NVSS has poorer resolution 
than FIRST, but greater sensitivity to extended structure. 
It also has a survey limit similar to FIRST, and of the 107 
SF galaxies in our sample with r > 10", 87 are present in 
the NVSS catalogue, within a 15" matching radius of the 
SDSS object (c.f. the 10" radius used by Sadler et al. 2002, 
for NVSS and 2dFGRS matches). These objects are shown 
in Figure 6, where the NVSS 1.4 GHz luminosities are com- 
pared with those from FIRST. For the systems larger than 
about 10" the extent of the underestimate in the FIRST 
measurements is clear. For all subsequent analysis herein, 
we use the NVSS derived SFRs in place of those from 
FIRST for galaxies with r > 10" to ensure no underesti- 
mates of 1.4 GHz derived SFR bias our results. Omitting 
the 20 objects with r > 10" and no NVSS measurement 
leaves a final sample of 771 SF galaxies. 
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3.2. Ha Star formation rates 

In estimating SFR from the Ha luminosity we adopt the 
calibration given by Kennicutt (1998): 

SFR Ha (M yO= i27 ^ o34w . (3) 

Prior to applying this calibration, though, there are several 
issues to address regarding the measurement of an Ha lu- 
minosity representative of the full emission from a galaxy 
given a flux measurement from a fiber-based spectrum. 
These issues include both corrections for obscuration, due 
to the intrinsic dust content of the target galaxy, and aper- 
ture corrections, to account for the emission missed by 
virtue of the fiber diameter potentially being smaller than 
a target galaxy. Both these effects can be accounted for 
using the available SDSS data. 

Before addressing these issues, we briefly digress to 
comment on the observing strategy called "smearing" 
(Stoughton et al. 2002). This is used in the SDSS to ac- 
count for seeing and wavelength-dependent atmospheric 
refraction effects in ensuring the most accurate spec- 
trophotometric calibration. A "smear-correction," in brief, 
consists of scaling a primary spectrum so its smoothed con- 
tinuum equals the smoothed continuum of a short expo- 
sure of the same object taken while the telescope pointing 
was dithered to cover a somewhat larger aperture than 
the fiber diameter. This scaling preserves line equivalent 
widths, but not line fluxes. The line fluxes are effectively 
being "aperture corrected" by scaling them assuming the 
line emission scales directly with the stellar continuum. 
This actually follows the same principle as the methods 
used below to perform aperture corrections, although the 
scaling of the smear-correction is to the stellar continuum 
within a fixed aperture for all galaxies, and will not neces- 
sarily encompass all of a target galaxy, especially at lower 
redshifts. 

Naively using measurements of spectral line fluxes from 
smear-corrected spectra can thus potentially lead to small 
errors in the resulting fluxes and line ratios. Because of the 
fixed aperture effectively used in constructing the smear- 
correction, the line fluxes from smear-corrected spectra 
will be different from those constructed from the primary 
spectra with an aperture correction based on measured 
galaxy sizes or magnitudes. Some preliminary compar- 
isons for line ratios, using Balmer decrement measure- 
ments in SF galaxies, quantify this small effect to some 
extent. Even with the Balmer decrement line ratio (widely 
spaced in wavelength) the systematic difference between 
primary and smear-corrected spectra is only about 10% at 
most. 

In summary, equivalent widths measured from smear- 
corrected spectra should be consistent with those in pri- 
mary spectra, and line ratios for lines close to each other in 
wavelength (such as the typical AGN/SF diagnostics) will 
be negligibly affected. But it is clear that to avoid poten- 
tially introducing small systematic offsets, the best solu- 
tion for individual flux measurements, or line ratios widely 
separated in wavelength, is to use spectral line measure- 
ments made on spectra which have been processed without 
using the smear correction. All spectral line measurements 
used for the analysis herein were made on primary spectra, 
with no smear corrections applied. 



3.2.1. Obscuration corrections 

Galactic foreground obscuration is corrected for follow- 
ing Schlcgcl et al. (1998), but obscuration by dust intrin- 
sic to the SF galaxies can cause more significant under- 
estimates in the emission line and u-band derived SFRs. 
We address this by making obscuration corrections in two 
ways. All the SF classified galaxies have measured Ha and 
H/3 fluxes, so the Balmer decrement Fna/Fn/3 can be cal- 
culated, and used to estimate and correct for the obscura- 
tion. The suggestion of luminosity-dependent obscuration 
in SF galaxies (e.g., Hopkins et al. 2001; Sullivan et al. 
2001), also recently identified in individual regions of star 
formation from spatially resolved spectroscopy of an ex- 
treme SF galaxy, IRAS 19254-7245 (Berta et al. 2003), is 
explored as well. New results from the Phoenix Deep Sur- 
vey (Afonso et al. 2003) indicate, however, that the situa- 
tion is more complex than implied by simple linear models. 
While galaxies with low SFRs seem to have relatively low 
levels of obscuration, at higher SFRs a broad range of ob- 
scurations are seen. A trend for an increase in the median 
obscuration with SFR, though, is still present. We make 
use of the method described by Afonso et al. (2003) and 
the relationship they derive to explore the effectiveness of 
such an empirical correction for the present sample, and 
perform an explicit comparison with the Balmer decrement 
correction. 

Stellar absorption in the Balmer emission lines, if not 
accounted for, can cause a significant overestimation in 
the implied obscuration from measurements of the Balmer 
decrement (Rosa-Gonzalez et al. 2002). Detailed analysis 
of continuum fitting to refine emission line measurements 
is being addressed by Tremonti et al. (2003), and an al- 
ternative line measurement method explored by Goto et 
al. (2003), but for the purposes of the current investiga- 
tion it was deemed sufficient to assume a simple constant 
correction for stellar absorption in the measured Balmer 
line equivalent widths (EWs). The value of this equivalent 
width correction, EW C , for the H/3 line was found by Miller 
& Owen (2002) to vary from about 1 A for Sa galaxies to 
about 4 A for extreme late types, with EW C = 2 A for Sb 
galaxies, consistent with the typical value found in other 
studies of SF galaxies (e.g., Georgakakis et al. 1999; Tresse 
et al. 1996). 

Given the broad range of possible stellar absorption val- 
ues, the assumption of a common value may act to intro- 
duce a degree of scatter into the resulting SFR estimates. 
For the sample investigated here, the stellar absorption 
correction is typically a significant fraction of the measured 
H/3 equivalent width, anything from a 10% to 100% correc- 
tion. As a result, it is worth emphasising that the extent 
of the derived obscuration correction depends strongly on 
the assumed value of the stellar absorption. The assump- 
tion of a common value as done here should therefore be 
restricted only to studies of large samples where the gross 
characteristics of the population are being examined, and 
more refined measurements should be preferred for analy- 
ses of individual objects. Furthermore, since the relative 
correction is largest for systems with low EW, and since 
these tend to be those with higher luminosities, it is the 
bright galaxies that are most affected by any uncertainty 
in the extent of the absorption correction. 
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The EW correction is converted into a flux correction 
for the Balmer lines using 



_ EW + EW C 
S - EW F ' 



(4) 



where S is the stellar absorption corrected line flux for 
Ha or H/3, EW is the equivalent width of the line, EW C 
the correction for stellar absorption, and F the observed 
line flux. The commonly assumed stellar absorption of 
EW C = 2 A at H/3 derives from a line measurement tech- 
nique which involves an integration of the spectrum over 
a wavelength range encompassing the line (Miller & Owen 
2002; Goto et al. 2003). When alternative measurements 
are used, such as the Gaussian fitting of the SDSS pipeline, 
it appears that a smaller stellar absorption correction is 
necessary. This is a result of the SDSS spectral resolution, 
which partially resolves the stellar absorption, resulting in 
the fitted Gaussian being only incompletely diminished by 
the absorption. Thus the measured line flux (or EW) is 
larger than a flux-summing method would produce. 

To establish the extent of this effect, the approximate 
distribution of actual stellar absorptions was first explored 
to determine the median absorption at the wavelength of 
Ha. Rather than simply assuming EW C = 2 A for H/3, 
we use the measured absorption at US in the SDSS spec- 
tra as a proxy for the H/3 absorption. We then followed 
the prescription of Keel (1983) for deriving the stellar ab- 
sorption at Ha, EW(Ha) = 1.30 + 0.40 EW(H/3) (Miller & 
Owen 2002). The resulting distribution of stellar absorp- 
tion EWs at Ha for a complete, volume limited sample 
of SF galaxies (constructed below, see §3.5) is shown in 
Figure 7. The median value of this distribution is 2.6 A, 
consistent with the observed range for late type galaxies. 
Since the resolution of SDSS spectra are sufficient to re- 
solve this absorption, a correction smaller than this actual 
value needs to be applied. This is demonstrated explicitly 
for an actual SDSS spectrum in the inset of Figure 7. The 
SDSS pipeline Gaussian fit to the emission line underes- 
timates the line flux by the black shaded area minus the 
grey shaded area, which can be seen to be smaller than the 
total stellar absorption, by at least a factor of two. Given 
the median stellar absorption of 2.6 A, it seems that the 
appropriate stellar absorption correction should be at most 
about EW C = 1.3 A. We adopt this value for the rest of 
this analysis. 

Now, using the stellar absorption corrected Balmer line 
measurements, the Balmer decrement Sn a /Sn/3 can be 
constructed. The obscuration correction is then derived 
using the Balmer decrement and an obscuration curve. For 
making obscuration corrections to the emission lines, we 
use the Milky Way obscuration curve of Cardelli et al. 
(1989) as referenced in Table 2 of Calzetti (2001). For the 
obscuration corrections to the stellar continuum in the it- 
band (§3.5), we use the obscuration curve of Calzetti et 
al. (2000) derived for starburst galaxies. Note that the ex- 
tent of the obscuration experienced by emission lines and 
the stellar continuum at the same wavelength differs by a 
factor of about two (Calzetti 2001). 

The stellar absorption corrected Balmer decrements are 
shown as a function of SFRi.4GHz in Figure 8. The pre- 
dicted Balmer decrement from the SFR-dcpendent obscu- 
ration of Afonso et al. (2003) is shown as the solid line in 
this Figure, and that of Hopkins et al. (2001) as the dashed 



line. The relation independently derived by Sullivan et al. 
(2001) lies 10% higher than that of Hopkins et al. (2001), 
while having an almost identical slope. These relations 
have been converted to the current SFR calibrations and 
cosmologies where necessary. As discussed by Afonso et 
al. (2003), the empirical relation found by Hopkins et al. 
(2001) appears to be affected by sample selection effects 
(being restricted to higher EW systems) . This results in a 
model that provides a reasonable approximation for larger 
EW systems (EW(Ha)> 70 A), but is clearly an underes- 
timate for smaller EW systems. For the whole sample, the 
model of Afonso et al. (2003) is better at tracing the trend 
in the median obscuration with SFR, although appears to 
be a slight overestimate. This is discussed further in § 3.2.3 
below. Again, very importantly, there is significant scat- 
ter in the distribution of actual Balmer decrements about 
such trends, particularly at high SFR. 

On top of the overall trend for systems with higher 
SFRs to have a higher median Balmer decrement, the 
observed trend with EW in Figure 8 is also likely to 
be real, and not merely an artifact of the chosen stellar 
absorption correction. In the relatively low SFR range 
1 < SFRi.4ghz < lOMoyr -1 , for example, the systems 
with EW< 70 A have a median Balmer decrement of 6.1, 
compared with 5.3 for the systems with larger EWs. If 
the stellar absorption correction was boosted as high as 
EW C = 2.6 A, (an unreasonably large estimate given the 
illustration in Figure 7), these median values decrease only 
to 5.4 and 4.6 respectively. This indicates both that there 
are real differences in the extent of the typical absorption, 
depending on the observed EW, and that this difference is 
not an artifact of the chosen stellar absorption correction. 
The extent of the absorption, moreover, is not negligible, 
even at these relatively low SFRs. 

3.2.2. Aperture corrections 

In addition to the obscuration correction, the emission 
line luminosities also require an aperture correction to ac- 
count for the fact that only a limited amount of emission 
from a galaxy is detected through the 3" diameter fiber. 
For both the Ha and [Oil] emission lines, this is done as 
follows. 

The Ha EW (corrected for stellar absorption) can be 
used along with an estimate of the continuum luminosity 
for the galaxy from the photometric catalogue at the ob- 
served wavelength of Ha, to recover an effective Ha line 
luminosity for the whole galaxy. Explicitly, (before obscu- 
ration corrections are applied), 

Lhq (W) = (EW + EW C ) l0-°- 4 (^- 34 - 10 ) ' ' 1 ' 



(6564.61(1 + z)) 2 
(5) 

where M r is the k-corrected absolute r-band AB- 
magnitude, derived from the observed r-band Petrosian 
magnitude. The last term converts this luminosity from 
units of WHz -1 to WA" 1 . This Equation assumes that 
the flux of the continuum at the wavelength of Ha can be 
represented by the flux at the effective wavelength of the 
r-band filter (« 6222 A). While this is not strictly true, the 
continuum in these SF systems in the wavelength range of 
interest is flat enough that it is a good approximation. A 
more refined estimate can be made using an interpolation 
between the absolute magnitudes in the r and i, or i and z 
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filters, as appropriate for the redshift of the galaxy. This 
was explored and found to make a negligible difference in 
the resulting distribution of SFRs and in the comparison 
with SFRs estimated at other wavelengths. 

The aperture correction implicitly assumes that the 
emission measured through the fiber is characteristic of the 
whole galaxy, and that the SF is uniformly distributed over 
the galaxy. Issues such as galaxy orientation and patchy 
distributions of SF regions, primarily when the aperture 
corrections are large, will increase the uncertainties in this 
form of aperture correction estimate. 

In the case of the [On] luminosity, an exactly analo- 
gous method is used to calculate the aperture correction, 
although no stellar absorption correction is necessary for 
the EW of [Oil]. In Equation 5 M u is substituted for M r , 
and the wavelength of [On], 3728.30 A, is used in place 
of that of Ha, 6564.61 A, in the last term. The explicit 
calculation used in constructing the final SFRh q > incorpo- 
rating both the aperture correction and the obscuration 
correction, is given as Equation B2 in Appendix B. 

An alternative method of applying an aperture correc- 
tion is described in Appendix A, along with a fuller dis- 
cussion of systematic uncertainties and an exploration of 
how the aperture correction depends on parameters such 
as angular size and redshift. 

3.2.3. Comparison with 1.4 GHz SFRs 

The estimated SFRs derived from Ha luminosities are 
shown as a function of SFRi. 4 gh z in Figure 9. Each step 
in the process of calculating the Ha luminosities is shown, 
to emphasise the magnitude of each effect. The uncor- 
rected SFRs shown are calculated directly from the Ha 
line luminosities prior to any obscuration or aperture cor- 
rections. The aperture corrected SFRs are calculated from 
the luminosities obtained using Equation 5 before apply- 
ing obscuration corrections, and finally, from Equation B2 
incorporating the obscuration correction as well. It can be 
seen that, on average, the combined effect of the aperture 
and obscuration correction is to increase uncorrected lumi- 
nosities (and SFRs) by a factor of about 20. The resulting 
estimates from Ha and 1.4 GHz luminosities are highly 
consistent, although there remains a significant amount of 
scatter. The rms characterising the extent of this scatter, 
in the sense of the rms deviation from the one-to-one line, 
is 0.21 dex, or about a factor of 1.6 either side of the line. 
The extrema of the scatter (apart from a small number of 
significant outliers) span about 1 dex. 

The results of using the obscuration correction method 
of Afonso et al. (2003) are also shown in Figure 9, produc- 
ing a similar distribution for the derived Ha SFRs as the 
Balmer decrement correction. Directly comparing SFRh q 
estimated using obscuration corrections derived from the 
Balmer decrement with those using the method of Afonso 
et al. (2003) show that the relation of Afonso et al. (2003) 
gives, on average, an overestimate of about 20%. This 
is not unexpected given the location of the trend shown 
in Figure 8. Afonso et al. (2003) discusses the fact that 
this relation, derived from a radio-selected sample of SF 
galaxies, reflects the presence of galaxies with larger ob- 
scurations, able to be detected in radio-selected samples, 
and possibly overlooked in UV or optically selected sam- 
ples. This point will be discussed in more detail in §4.2 



below. In any case, for this type of application empirical 
relations like that of Afonso et al. (2003) and Hopkins et 
al. (2001), when their selection effects and limitations are 
correctly taken into account, appear to be useful tools for 
estimating obscuration corrections in the absence of more 
physical measurements. 

It is possible that some of the scatter observed in Fig- 
ure 9 may be produced by our assumption of a common 
stellar absorption correction. This was investigated by ap- 
plying corrections using the individual estimates derived 
above, from the method of Keel (1983). No measurable 
reduction in the resulting scatter was detected, suggesting 
that the assumption of a common EW correction for stel- 
lar absorption does not dominate the observed scatter in 
these results. To account for the scatter, other physical 
processes must be investigated. 

The differences between the measurements of SFR from 
Ha and 1.4 GHz luminosities from Figure 9(c) were fur- 
ther investigated to explore whether there was any residual 
trend in the relation. Figure 10 shows the ratios of these 
SFRs as a function of redshift, indicating the consistency 
of the two measurements, on average. The effect of the flux 
density limit of the radio survey is also shown. This limit 
(about 0.75 mJy, somewhat lower than the completeness 
level of 2 mJy) starts to bias the observed distribution for 
systems beyond a redshift of z 0.1. This bias is in the 
sense of losing sensitivity to systems with apparent under- 
estimates of SFRi.4GHz with respect to SFRh q . Also, be- 
low a redshift of z m 0.03, the distribution seems to show 
predominantly more systems with SFRr- q overestimated 
compared to SFRi.4GHz- This is unlikely to be a selection 
effect, and is also unlikely to be an underestimate in the 
1.4 GHz luminosities, as the NVSS measurements used for 
these large, nearby systems are sensitive to the emission 
over the full range of sizes seen. This effect is instead more 
likely to be an overestimate in the Ha luminosities as a re- 
sult of the effective aperture correction used. The aperture 
correction implicitly assumes that the star formation sam- 
pled through the fiber is representative of the distribution 
over the whole galaxy, and that this scales directly with the 
broadband optical emission. The more realistic scenario is 
that (especially in relatively low SFR, nearby systems) the 
star formation is patchy, and distributed non-uniformly 
throughout the galaxy. Hence it is not surprising that the 
simple aperture correction used for these systems produces 
an overestimate. This overestimate is not seen in the dis- 
tribution of measurements between 0.05 < z < 0.1, where 
the 1.4 GHz flux density limit does not yet affect the ob- 
served distribution. The most likely explanation here is 
that the lowest redshift systems have a larger contribution 
from galaxies with centrally concentrated star formation, 
and possibly also from irregular and dwarf systems, known 
to contain patchy SF distributions. 

It is worth pointing out here that there is a subtle selec- 
tion effect arising from the requirement of having a radio 
detection for the SF sample. This is an effective emission 
line flux limit resulting from the radio flux density limit. 
If the 1.4 GHz detection limit of about 0.75 mJy is trans- 
formed to an Ha flux limit via the SFR calibrations of 
Equations 1 and 3, (after incorporating the factor of 20 
corresponding to the average combined obscuration and 
aperture correction), an effective Ha flux limit of about 
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3.5 x 10~ 18 Wm -2 is derived. This is over an order of 
magnitude higher than the minimum observed Ha fluxes 
of w 3 x 10~ 19 Wm~ 2 for the spectroscopically classified 
SF galaxies in the whole DM. A similar relation can be 
derived for [On] line fluxes. In this case the effective [On] 
flux limit («7x 10~ 19 Wm~ 2 ) imposed by the radio flux 
density limit is a little closer to the actual minimum [Oil] 
fluxes observed (wlx 1 0~ 19 Wm~ 2 ) in the DR1 SF galax- 
ies. The primary result of this selection effect seems to be 
that the radio detected sample is more similar to a com- 
plete sample, in the sense that the incompleteness at low 
Ha fluxes is removed. This is explored further in §4.2. 

3.3. [On] Star Formation Rates 

An SFR calibration based on the [On] emission is an 
important tool for probing galaxy SFRs at z > 0.4 where 
Ha is rcdshiftcd out of the bands easily accessible by op- 
tical spectroscopy. As with Ha, [On] luminosities require 
aperture and obscuration corrections, which are described 
in detail in § 3.2 above. SFRs derived from [On] lumi- 
nosities are based on the fact that there is a good cor- 
relation between observed [On] line fluxes and observed 
Ha lines fluxes (i.e., prior to any obscuration corrections). 
The SFRs come from using the [On] emission as a proxy 
for the Ha emission, which must then be corrected for ob- 
scuration based on the obscuration at the wavelength of 
Ha in order to use the SFR calibration derived for Ha 
luminosities (Kennicutt 1998, 1992). 

Many recent estimates of [On] SFRs rely on the SFR cal- 
ibration of Kennicutt (1998), an average of the calibrations 
reported by Gallagher et al. (1989) and Kennicutt (1992), 
each of which are based on samples of fewer than 100 galax- 
ies. In particular the observed ratio FiQu]/F Ka = 0.45 
from Kennicutt (1992) is extensively used. This ratio is, 
however, luminosity and metallicity dependent (Jansen et 
al. 2001), and the appropriate ratio should be determined 
for a given sample depending on its selection effects. The 
current sample of 791 SDSS SF galaxies with radio coun- 
terparts includes 752 with measured [On] emission. Fol- 
lowing the methodology used in Kennicutt (1992), the ra- 
tio of the observed fluxes in the [On] and Ha emission 
lines was investigated (see also Figure 20, and subsequent 
discussion in §4.2), and the median ratio was found to 
be F[ou]/Fu a — 0.23. (This result comes after a stellar 
absorption correction EW C = 1.3 A is applied to the Ha 
emission, although the median ratio changes negligibly if 
this step is omitted.) 

Applying or omitting the aperture corrections when cal- 
culating this ratio also changes the median value negligi- 
bly, although it is certainly true that a differential dis- 
tribution for the origin of the line emission would not be 
accounted for with the aperture correction methods used 
here. If [On] emission came predominantly from galaxy 
disks, and Ha predominantly from galaxy nuclei, for ex- 
ample, the current estimate would be skewed low. It seems 
unlikely, though, that this contrived geometry of line emis- 
sion should occur for the majority of systems. Given, addi- 
tionally, the expected association of both forms of emission 
with star formation regions, it seems reasonable that both 
[On] and Ha emission should be at least approximately 
colocated throughout SF galaxies. 

The median ratio determined for the current sample 



(from fiber spectroscopy) is slightly lower than that found 
by Jansen et al. (2001) (using spatially integrated spectra) 
for galaxies of similar absolute magnitude. This does not 
appear to be an artifact of the radio selected nature of the 
current sample. If the radio detection requirement is re- 
laxed, and a complete, volume-limited sample of optically 
selected SF galaxies constructed, an unbiased estimate of 
the distribution of this flux ratio can be established. This 
is pursued further in § 4.2 below, and results in median val- 
ues for Fjonj/Ffja very similar to that found for the full 
radio detected sample currently being explored. For the 
present discussion we adopt the ratio i^on] /-Fka = 0.23 
as being representative of typical SF galaxies, given the 
optical luminosity range of the present sample. This cal- 
ibration may also be useful for higher rcdshift surveys of 
galaxies of similar luminosity, where observations of the 
Ha line are more difficult. The metallicity dependence 
of the F[oii]/-fHa ratio, however, must also be accounted 
for at higher redshifts, given the strong evolution in the 
metallicity-luminosity relation (Kobulnicky et al. 2003). 

Using this estimate for the correlation between [On] and 
Ha line fluxes the calibration of Lrq to SFR (Equation 3) 
is transformed to 

SFR [on] (M yr-^^2L, (6) 

where £[oii]> due to the way this calibration was de- 
rived, must incorporate the obscuration correction valid 
at the wavelength of Ha (Kennicutt 1998). The final [On] 
SFR estimate is explicitly shown in Appendix B as Equa- 
tion B5. 

The effectiveness of this calibration can be seen in the 
comparison of [On] derived SFRs with those from 1.4 GHz 
and Ha in Figure 11. The small number of systems in 
Figure 11(a) with anomalously high 1.4 GHz SFRs may be 
composite AGN/SF systems, in which the AGN is masked 
by obscuration at optical wavelengths, since the Ha and 
[On] SFRs for these systems are in good agreement. 

3.4. FIR Star Formation Rates 

FIR flux densities are available from the IRAS catalogs 
for a sub-sample of our SF galaxies, and positional match- 
ing with the FIRST radio sources identifies 191 galax- 
ies in this final SF sample with IRAS detections. The 
60 fim and 100 /im flux densities, Sqo and Sioo can be 
used to derive a FIR flux following Hclou et al. (1988), 
FIR(Wm- 2 )= 1.26 x 10~ 14 (2.58 S 60 + S wo ). There are 
several calibrations of FIR luminosity to SFR in the lit- 
erature, (for a summary see Kennicutt 1998). We begin 
with the recent calibration of Bell (2003), (from which our 
chosen 1.4 GHz SFR calibration was derived): 

SFRfir (Mq yr ) - 2 44 x lo36 w " L39 x 10 36 W - 

(7) 

where 

f= ( 1 + y/2.186 x 10 35 W/£fir £fir > L c 

I 0.75 (1 + ^2.186 x 10 35 W/£ F ir) L F ir < L c , 

and L c = 2.186 x 10 W. Lfir is the luminosity corre- 
sponding to the FIR flux as defined above, and the factor 
of 1.75 in Equation 7 converts this to a luminosity repre- 
sentative of the full (8 — 1000 //m) mid- to far- infrared spec- 
trum (for details see Bell 2003; Kcwley et al. 2002). The 
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piecemeal nature of this Equation comes from the assump- 
tion of varying fractions of old stellar populations with lu- 
minosity. It is plausible, however, that galaxies with high 
FIR luminosities in the current sample have similar rel- 
ative contributions (to FIR luminosity) from old stellar 
populations as the lower luminosity systems. This is in 
contrast to the sample from Bell (2003), which was more 
inhomogeneously constructed, possibly resulting in prefer- 
ential selection of starbursting galaxies with low old frac- 
tions above 10 11 M Q yr _1 (Bell, 2003, private communica- 
tion) . Consequently we have chosen to adopt an old stellar 
population fraction of 30% independent of luminosity, and 
apply the portion of Equation 8 valid for Lfir < L c to all 
galaxies. This is explicitly presented in Equation B10. We 
further recommend this strategy to others working on sam- 
ples of radio-selected galaxies (Bell, 2003, private commu- 
nication). The resulting SFRfir is compared with those 
from 1.4 GHz and Ha in Figure 12. This calibration is 
within a factor of two of that of Kennicutt (1998) down to 
about 3 x 10~ 3 L c (Bell 2003). Note that the empirical re- 
lation from Buat & Xu (1996), measured from a sample of 
galaxies of types Sb and later, produces SFRs almost 80% 
larger than this calibration, and that of Condon (1992) 
produces SFRs about 60% larger. 

3.5. u-Band Star Formation Rates 

While ultraviolet (UV) luminosity (A < 2500 A) is com- 
monly used as an SFR indicator, the luminosity at u-band 
wavelengths (A s=s 3600 A) is similarly dominated in star- 
burst galaxies by young stellar populations, and in the 
absence of UV measurements the w-band luminosity may 
thus be used instead as an SFR indicator (Cram et al. 
1998). The u-band luminosities used here are derived from 
the SDSS k-corrected absolute u-band magnitudes (Blan- 
ton et al. 2003), after incorporating an obscuration cor- 
rection based on the Balmer decrement and the extinction 
curve from Calzctti et al. (2000). It is worthwhile pointing 
out that the average w-band obscuration correction ranges 
from a factor of 3 at SFRs of lM yr _1 up to about a 
factor of 10 at SFRs of lOOMoyr- 1 . 

UV luminosities have been extensively used as SFR in- 
dicators, and for wavelengths 1500 A < A < 2500 A the 
calibration given by Kennicutt (1998) 

has proven quite effective. For w-band luminosities, 
though, it is somewhat more difficult to assign a simple 
scaling factor to derive an SFR due to the strong de- 
pendence on the evolutionary timescale. From synthetic 
galaxy spectra (e.g., Fioc & Rocca-Volmerange 1997) it 
can be seen that the w-band luminosity varies from about 
factor of 10 lower than the UV luminosity at the onset 
of a burst of star formation, to almost equivalent by 10 8 
years later (see also discussions of the dependence of UV 
measures on the timescale of SFR in Sullivan et al. 2000; 
Glazcbrook et al. 1999). 

Given this sensitivity of the u-band (and indeed UV) 
luminosity to the starburst age and the assumed star for- 
mation history, a more complex calibration is in general 
likely to be necessary. This may take the form of a non- 
linear dependency on L\j to reflect the rapid change of 
L\j with respect to iuv during the first 10 8 years of a 



starburst, and to account for the presence of old stellar 
populations that are likely to contribute significantly to 
L\j in less luminous systems (Bell 2003). More quiescent 
or low-luminosity SF systems have a relatively larger con- 
tribution to their w-band luminosity from old stellar popu- 
lations, (they are, on average, redder than more luminous 
galaxies), and this causes a linear calibration from lumi- 
nosity to SFR to result in an overestimate of the SFR. 
While these effects can be modelled in detail using stellar 
spectral synthesis methods (Sullivan et al. 2000), for the 
purposes of a simple SFR estimate based on a measured 
luminosity an empirical non-unity power-law relationship 
between u-band luminosity and SFR can be derived. Al- 
though providing no information about the SFR history or 
the SED evolution, such a calibration has the advantage 
of being easy to construct and can form a useful tool in 
subsequent analysis. 

To construct such a tool we require a complete, volume- 
limited sample of galaxies, a necessity not demanded of 
the previous sections where straightforward comparisons 
between existing calibrations were being explored. This 
requirement is necessary to avoid biasing the calibration 
low or high based on the independent limits of the Ha 
and u-band detections. In order to have the largest com- 
plete sample from which to derive the new calibration, 
and to eliminate any potential concerns introduced re- 
garding radio selection, we thus construct an SF galaxy 
sample from the DR1 based on the completeness crite- 
ria described by Gomez et al. (2003), and do not re- 
quire radio detection at all. The criteria used specify that 
galaxies should lie closer than z < 0.095 (since aperture 
corrections are applied, no lower rcdshift limit is neces- 
sary) and have M T < —20.57 (after converting to our 
chosen value of H ). The small number of systems with 
z < 0.05 (the lower redshift limit used by Gomez et al. 
2003) which show slight overestimates in their aperture 
corrected SFRna relative to SFRi.4qhz do not bias the 
result, which remains unchanged if these objects are ex- 
cluded. A further restriction was applied, a limit on the 
Ha fluxes, F Ka > 3 x 10~ 18 Wm~ 2 . This corresponds 
to an SFR of about 0.5M Q yr _1 at the redshift limit of 
the complete sample, and ensures no bias will be intro- 
duced to the derived calibration through the presence of 
incompletely sampled fainter Ha systems. This results in 
a sample of 2625 spectroscopically classified SF galaxies 
with measured Ha and u-band luminosities. 

Using this sample the ordinary least-squares bisector 
method of linear regression (Isobe et al. 1990) was ap- 
plied to log(Lu) (after applying obscuration corrections 
based on the Balmer decrements) and log(SFRn Q ). It is 
recognised that the luminosity limits imposed by requiring 
completeness will cause a small bias to this fit, (since the 
imposed limit on r-band luminosity translates to an effec- 
tive limit on w-band luminosity, and omitting galaxies with 
lower luminosities affects the slope of the fit slightly), nev- 
ertheless this still produces a better estimate of the true 
relation than were an incomplete sample used. The com- 
parison for the complete sample of galaxies between L\j 
and SFRh q can be seen in Figure 13(a), which indicates 
the resulting best fit relation. This relation results in the 
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calibration: 

/ T x 1.186 

^^^ = ( iCTw^ J ' (10) 

and the result of applying this calibration is shown in Fig- 
ure 13(b). Strictly speaking, this calibration is only valid 
over the range of luminosities probed here, 2 x 10 21 < 
Ljj < 10 23 WHz _1 . The explicit calculation in terms of 
the observable quantities is given in Equation B8, in Ap- 
pendix B. Since the calibration from luminosity to SFR 
is no longer linear, a multiplicative obscuration correction 
cannot be applied directly to the SFR in the same manner 
as to the luminosity. Rather the appropriate exponent as 
given in the SFR calibration must be incorporated as well. 

The entire sample of 21649 SF galaxies with measured 
u-band and Ha SFRs are compared in Figure 14. The 
effects of the incompleteness can be seen, appearing as a 
slight apparent bias to overestimated w-band SFRs. The 
fitted relation, therefore, is sensitive to the completeness of 
the sample being used, emphasising the need for the initial 
calibration to be performed on a well-defined sample. 

Using the new SFR calibration for u-band luminosity, 
SFRu is shown as a function of SFRi.4GHz in Figure 15. 
The rms deviation between SFRu and SFRi. 4 gh z is 0.23 
dex, or a factor of 1.7 either side of the one-to-one line. 

4. DISCUSSION 

4.1. Systematics and absolute SFR calibrations 

We have demonstrated a consistency between indepen- 
dent SFRs estimated from Ha and 1.4 GHz luminosities, 
used an observed Ha/ [On] flux correlation to establish the 
consistency of the [On] SFRs, identified the FIR calibra- 
tion giving most consistent SFRs, and derived a u-band 
calibration defined to give consistent SFRs. It should be 
emphasised that, since the calibration for SFRi.4GHz is 
derived from the FIR calibration (Bell 2003), and since 
the [On] and u-band calibrations are derived from the 
Ha calibration, the only independent calibrations being 
directly compared are those of Ha and FIR, both ulti- 
mately coming from Kennicutt (1998). This is still non- 
trivial, however, given the different physics involved in the 
two derivations (ionising luminosities compared with total 
bolometric fluxes), and the broad consistency seen across 
the range of wavelengths explored is highly encouraging. 
This consistency is in spite of the large scatter, frequently 
referred to above, in the SFR values measured at different 
wavelengths for individual systems. The dispersions about 
the one-to-one line for each estimator compared with the 
1.4 GHz estimator are summarised in Table 1, where the 
uncertainties shown are the rms deviations either side of 
the one-to-one line. Some issues that are likely to con- 
tribute to this observed scatter include (1) the reliability 
of the assumptions underlying the aperture corrections to 
the emission line estimates; (2) whether the Balmer decre- 
ment obscuration estimate (which is measured only for 
the region seen through the fiber aperture) is valid for the 
whole galaxy; (3) the extent to which any low luminosity 
AGN that may be present in the SF galaxies contributes 
to the observed radio luminosity; (4) galaxy to galaxy dif- 
ferences in the average 1.4 GHz luminosity generated by 
supernovae, either through differences in electron popula- 
tion densities, temperatures, or confinement, strength of 
interstellar magnetic fields, or other physical differences. 



It is possible that in addition to this observed scat- 
ter, however, the absolute SFR calibrations are still un- 
certain by up to a factor of two (Bell 2003; Condon 2002). 
Apart from concerns arising through the detailed physics 
involved in deriving the SFR calibrations, the only ad- 
ditional systematics that might change the quantitative 
SFRs derived herein are the assumed equivalent width cor- 
rection to account for stellar absorption, and the chosen 
obscuration curve. The aperture corrections could produce 
underestimates of SFRh q only for galaxies large enough 
that the fiber predominantly samples light from the nu- 
cleus, the star formation occurs primarily in the disk, and 
the galaxy is observed close to face on. Hence it seems 
likely that any bias caused by the aperture corrections will 
be overestimates of SFRh q , as seen in Figure 10, and that 
these are relatively small effects restricted to very nearby 
systems. (It is of course possible that the assumptions in- 
volved in making the aperture correction contribute, per- 
haps significantly, to the scatter seen in the measured in- 
dividual SFRs.) 

Differences between estimates of the obscuration curve 
appropriate for correcting emission lines produce only 
small changes (of order 5%) from the results given here. 
The assumption of Case B recombination may also con- 
tribute to the scatter between the different SFR estimates, 
since this will not be valid for all systems (although it 
should be reasonable for the majority of the SF galax- 
ies). It is also possible that the detailed geometries of 
gas and dust in individual objects could result in quite 
different obscurations than estimated here, but for this 
sample of relatively low redshift galaxies the obscuration 
curve adopted should be fairly representative. These ef- 
fects, while contributing to the observed scatter, should 
not act in a systematic fashion. 

The Ha SFR values reported here could potentially be 
increased by up to about 20% by reducing the chosen EW 
correction for stellar absorption (a value of EW C = 0.7 A 
increases the derived SFRs by this factor). This arises 
through the change introduced in the Balmer decrement 
when a different value of the EW C is used. Since the mea- 
sured H/3 emission is in many cases comparable to the cor- 
rection being applied to it, the resulting Balmer decrement 
is very sensitive to the chosen value. Values of EW C smaller 
than about 0.7 A would, however, no longer be consistent 
with the stellar absorption corrected Balmer decrements 
derived from the flux-summing method of line measure- 
ment. Making the stellar absorption correction smaller 
also increases the discrepancy between the two methods 
(described in Appendix A) of estimating SFRh q . So there 
is a limit of about 20% by which SFRjja can be increased 
through the chosen estimate of EW C . Even in combination 
with the other systematic uncertainties in the Ha SFR 
calibration, the values reported here seem reliable in an 
absolute sense to better than a factor of two. In any case, 
it doesn't seem to be possible to produce high enough val- 
ues of SFRjja to become consistent with the 1.4 GHz SFR 
calibration of Condon (1992) (which gives SFRs about a 
factor of two higher than the calibration used here), with- 
out revising the Ha SFR calibration. 

Since the SFR[qii] calculation uses the same obscura- 
tion correction as for Ha, and since the u-band SFRs are 
calibrated to SFRh«, these would both be changed con- 
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sistently as different estimates of EW C change SFRh q in 
the above case. The FIR estimate, similarly, is consis- 
tent with the chosen 1.4 GHz estimate (by construction), 
so the only truly independent SFR estimators in the cur- 
rent investigation are SFRh q and SFRi.4ghz- The con- 
sistency between these two results does, however, support 
the reliability of the assumed stellar absorption correction, 
and since the absolute calibration of the Ha SFRs appears 
quantitatively reliable, to better than a factor of two, the 
radio and FIR calibrations adopted must be similarly reli- 
able. This conclusion regarding the quantitative reliability 
of the SFR calibrations, of course, only applies to the over- 
all population. Measurements of individual objects, as can 
be seen from the large scatter in all the SFR comparison 
diagrams, can still easily be uncertain by factors of two or 
more. 

So, despite the discrepancies (factors of two) in many of 
the available calibrations, a strong argument can be made 
that the absolute SFR calibrations should be close to those 
adopted here. This preliminary exploration suggests that 
the presently adopted calibrations give SFRs which are on 
average reliable to better than a factor of two. 

4.2. Radio-selected SF galaxies 

The properties of the radio detected SF galaxies have 
been explored to investigate some details of how they differ 
from optically selected SF galaxies. The complete sample 
of SF galaxies constructed for the it-band calibration in 
§3.5, above, was used to compare the distribution of var- 
ious measured parameters between the optically selected 
and radio detected objects. There are 2625 galaxies in 
the complete sample, of which 380 are radio detected. 
Figure 16 indicates the distribution in optical luminos- 
ity spanned by both the complete sample and the radio 
detected galaxies within that sample. The median magni- 
tudes of these two distributions are similar (—20.6 for the 
complete sample, —20.8 for the radio detected systems), 
although the radio detected systems show a much broader 
magnitude distribution. Gaussian fits to the two his- 
tograms indicate a full width at half maximum (FWHM) 
of 0.9 mag for the complete sample, and 1.3 mag for the 
radio detected galaxies. Figure 17 shows the distributions 
of rest frame u — g and u — r colors, concentration index in- 
verse, 1/c, and D4000 for both the complete SF sample and 
the radio detected galaxies in the complete sample. (Con- 
centration index, here, is defined to be the ratio of the 
Pctrosian radius enclosing 90% of a galaxy's light to that 
enclosing 50%, e.g., Kauffmann et al. 2003.) It is clear 
that the radio detected galaxies preferentially identify a 
sub-population of the SF galaxies, with somewhat redder 
colors, more concentrated morphologies, and higher D4000 
values. This all suggests that the radio selected SF galax- 
ies favour systems with a greater relative contribution from 
old stellar populations than optically selected systems. It 
should be noted here that while the radio detected SF 
systems are redder and bulgier than optically selected SF 
galaxies, they are still on average bluer and diskier than 
the overall population of galaxies. This is confirmed by the 
dashed histograms in Figure 17, which show the distribu- 
tion of parameters for a complete sample of all galaxies, 
constructed in an identical fashion to the complete sample 
constructed above, but without any restriction on the Ha 



EW. This new complete sample of 24444 galaxies clearly 
shows bimodal distributions corresponding to red and blue 
galaxy populations, and verifies that the radio-detected SF 
galaxies truly represent a sub-population of SF galaxies, 
and are not the result of an unexpected selection effect, or 
some unexplained early-type galaxy contamination to the 
SF sample. 

The broader distribution seen in absolute magnitude for 
the radio detected galaxies suggests that these results may 
simply reflect a tendency for the radio detection to favour 
larger, brighter galaxies. But the distributions in the it- 
band luminosities, for the radio detected galaxies com- 
pared with the complete sample, suggests that this may 
not be the whole picture. In Figure 18 the u-band and 
Ha SFRs are compared for the galaxies in the complete 
sample, as for Figure 13(b), but now identifying the ra- 
dio detected galaxies. There seems to be a preference for 
the radio detected systems to appear below the one-to-one 
line, at least for low to moderate SFRs. Since the compar- 
ison between radio and Ha SFRs shows no such effect, it 
seems likely that the radio detection preferentially identi- 
fies lower u-band luminosity systems for a given SFR. 

The reasons why these types of preferential selection 
occur may be related to the optically selected samples un- 
dersampling the red end of the distribution as a result of 
obscuration effects, but it is also possible that the radio de- 
tection may undersample the blue end of the distribution 
due to the different physical processes and timcscales pro- 
ducing the emission in the different wavelength regimes. 
A full exploration of the questions raised by these and re- 
lated results is clearly warranted, although, being beyond 
the scope of the present work, this is being investigated in 
detail in a subsequent paper (Hopkins et al. 2003). 

The radio detected systems seem to have a notably 
higher median obscuration than the complete sample, 
ijia = 1.6 mag compared with An a = 1.2 mag (Fig- 
ure 19), and the distribution is broader, a Gaussian fit 
giving a FWHM of 1.2 mag compared with 1.0 mag for the 
optically selected galaxies. The median value here for the 
optically selected objects is consistent with the 1.1 mag of 
obscuration commonly assumed for Ha measurements of 
SF galaxies (Kennicutt 1983). The median value for the 
radio detected systems, however, is significantly higher, 
and may be related to the efficiency of radio measurements 
in detecting heavily obscured sources. In other words, 
a radio selected sample should not bias against galaxies 
with high obscuration, producing a distribution of obscu- 
rations characteristic of the true distribution (Afonso et al. 
2003). The slightly higher median and broader distribu- 
tion for the radio detected systems suggests that it may be 
likely that the optically selected sample does indeed under- 
sample the redder, bulgier end of the galaxy distribution. 
This docs not, however, exclude the possibility that ra- 
dio detection simultaneously undersamples the bluer end, 
as suggested by the lack of higher SFRu systems seen in 
Figure 18. 

In Figure 20(a), histograms of the Ffonj/^Ha flux ra- 
tio for the (incomplete) sample of all the SF galaxies from 
DR1 shows a median value of 0.38, similar to that found 
by Kennicutt (1992). The radio selected systems here 
have a much lower ratio, 0.23, more comparable to the 
median values found for the complete samples, shown in 
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Figure 20(b), for both optically selected and radio detected 
galaxies (0.26 and 0.21 respectively). It appears, moreover, 
that there is a weak trend between the ratio Fpuj/^Ha 
and Ha line flux or EW, the latter being shown in Fig- 
ure 21, in the sense of higher flux ratios in higher EW 
systems. This is likely to be directly related to the corre- 
lation with luminosity noted by Jansen et al. (2001), as the 
higher EW systems include those with the lower luminosi- 
ties. The discrepancy between the present result and those 
from Gallagher ct al. (1989) and Kennicutt (1992) can also 
be explained by some level of incompleteness in those sam- 
ples at the low EW end. This would lead to a bias toward 
higher EW systems, and hence higher F[oirj/^Ha ratios. 
For systems in the complete sample with EW(Ha)> 70 A 
the mean ratio is 0.46, much closer to the mean of 0.45 
found by Kennicutt (1992). 

An important conclusion here is that the use of [On] 
luminosities as an SFR estimator requires a good under- 
standing of the selection effects of the sample being in- 
vestigated. The notion that SFRs based on [On] lumi- 
nosities are not very reliable (Jansen et al. 2001; Char- 
lot & Longhctti 2001; Kennicutt 1998) is perhaps related 
at least partially to the incomplete understanding of the 
F[ou] /Fho distribution valid for the sample under consid- 
eration, and what governs this distribution for the given 
set of selection effects. To the extent that the Froii]/FH Q 
ratio depends on mctallicity, for example, the evolution 
of the luminosity-metallicity relation (Kobulnicky et al. 
2003) must also be taken into account for higher redshift 
systems. 

A brief exploration of how Fjoirj/Fka varies with several 
parameters was pursued, to identify any obvious trends 
that might help in this type of sample characterisation. 
Jansen et al. (2001) find a strong correlation between 
i^onj/^Ha and Mb, with brighter galaxies having lower 
flux ratios. No such relation (with M g , M r or M z ) is 
convincingly measurable in the present sample although 
this only spans a range of about three magnitudes, much 
smaller than the eight magnitudes spanned by the sample 
of Jansen ct al. (2001). Given the scatter in the relation 
shown in Figure 1 of Jansen et al. (2001), a range of at 
least 5 or 6 magnitudes would need to be sampled to iden- 
tify this trend, so it is not surprising that it is not clearly 
detected in the current sample. There is, however, a sug- 
gestion in the current data that the brightest systems are 
restricted to lower ratios, and this weak trend is reflected 
in the weak trend mentioned above with Ha EW. There 
is, similarly, very little detected trend with galaxy size, 
although again the very largest galaxies are restricted to 
lower ratios. 

The connection between the F[oii]/Fh q ratio and the 
B aimer decrement was also examined. If a common obscu- 
ration curve is valid for all the galaxies in the sample, there 
should be a simple relationship between this ratio and the 
Fkc/Fhp ratio. Figure 22 shows this relation for the com- 
plete sample of SF galaxies (compare with Figure 2(c) of 
Jansen et al. 2001). The curves shown in Figure 22 in- 
dicate the relationship expected for the obscuration curve 
of Cardelli et al. (1989) and intrinsic F] OII ]/Fh q flux ra- 
tios (before any obscuration affects the emission lines) of 
0.3, 1.0 and 2.0. This suggests either that the complete 
sample displays a range of these intrinsic ratios spanning 



about 0.3 — 2.0, or that the attenuation in SF galaxies does 
not follow a unique attenuation law owing to varying dust 
geometries, changes in dust properties, or both. 

Another comment can be made regarding Figures 21 and 
22. A higher value for F|oh]/Fh q corresponds to lower ob- 
scuration (from Figure 22), and EW(Hq) is a tracer of the 
current relative SF (in the sense of absolute current SFR 
relative to the total integrated SFR). So the weak trend 
seen in Figure 21 might be taken to suggest that systems 
with higher relative SFRs show (on average) lower obscu- 
rations. This is an intriguing result given that Figure 8 
implies systems with higher absolute SFRs have (on aver- 
age) higher obscurations. These results may be useful in 
further exploring the nature of obscuration in SF galaxies 
as a function of their physical properties. One hypothesis 
which will be examined in more detail in a subsequent in- 
vestigation is that these results are both consistent with 
obscuration being primarily a function of galaxy mass or 
luminosity (larger galaxies having higher obscurations). 

5. SUMMARY 

From a sample of 3079 SDSS galaxies having radio lu- 
minosities from the FIRST survey, we have used opti- 
cal spectroscopic diagnostics to identify a sub-sample of 
791 SF galaxies. Using this sub-sample we have investi- 
gated five SFR indicators based on Ha, [On], u-band, FIR 
and 1.4 GHz luminosities. The FIRST 1.4 GHz derived 
SFRs below about 10M Q yr _1 are progressively underes- 
timated as the galaxy sizes become larger, consistent with 
a known limitation of the FIRST survey. This can be 
corrected by using NVSS measurements where available 
for the larger systems. After applying appropriate obscu- 
ration and aperture corrections, the Ha SFR estimate is 
seen to be consistent with the 1.4 GHz estimate, although 
a large scatter still remains. The median [On] to Ha flux 
ratio is found to be F[ OII ]/Fh q = 0.23, about a factor of 
two lower than commonly assumed, and an updated SFR 
calibration for [Oil] luminosities was derived to account 
for this. The resulting [On] SFRs are highly consistent 
with those from Ha and 1.4 GHz luminosities. A power- 
law calibration between u-band luminosities and SFRh q 
was empirically derived, which provides measurements of 
SFR based on L\j consistent with the other estimators, 
although more detailed investigation of the physical pro- 
cesses driving this relation are warranted. Issues surround- 
ing the reliability of the absolute SFR calibration were ad- 
dressed, suggesting that the calibrations used here should 
be reliable to better than a factor of two. 

With these results there are now three reliable SFR esti- 
mators available from the SDSS measurements, Ha, [On], 
and u-band luminosities. With the large sample size and 
the extensive spectroscopic measurements archived by the 
survey, it thus provides a vast resource for investigations 
of star formation in the universe. 

Investigating the properties of the SF galaxies, it is 
found that the median obscuration at the wavelength of 
Ha for the complete sample of SF galaxies is Ah q = 
1.2 mag, comparable with the 1.1 mag from Kennicutt 
(1983), while the radio detected systems are notably 
higher, Ah q = 1.6 mag. The properties of the radio 
detected sample imply that radio selection preferentially 
identifies somewhat redder, bulgier SF systems, (although 
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still bluer and diskier than the general galaxy population), 
having a relatively larger contribution from the old stel- 
lar population, than seen in optically selected SF samples. 
This is characterised through redder optical colors, and 
higher D4000 values and concentration indices than in op- 
tically selected samples. This could be attributed to either 
or both of the cases that (1) the optically selected samples 
undersample the red end of the distribution due to obscu- 
ration effects, and (2) radio detection undersamples the 
blue end of the distribution due to the different physical 
processes and timescales producing the multiwavelength 
emission. 
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APPENDIX 

A. APERTURE CORRECTIONS 

In addition to the aperture correction described by Equation 5, the Ha luminosity calculated using the line flux can 
be explicitly aperture corrected. This is done using the ratio of the fluxes corresponding to the total galaxy magnitude, 
and the magnitude "through the fiber," where this latter term is also one of the outputs of the photometric pipeline. 
This fiber magnitude comes from a photometric measurement of the magnitude in an aperture the size of the fiber, and 
is corrected for seeing effects. The Petrosian magnitude was used to represent the total galaxy flux. The extent of this 
aperture correction can be expressed as 

A = 10- - 4 (^tro-r fib er) ( Al ) 

where r^bor is the r-band fiber magnitude. The explicit calculation of the aperture corrected luminosity using this method 
is thus 

Ln a (W) = AnDf SW 10- a4 (^ tr o-r flbor ) _ ^ 

with Sn a being the stellar absorption corrected flux of the Ha emission, and Di the luminosity distance. The obscuration 
correction has not been included in the above Equation. 

Both methods for estimating an aperture-corrected emission line luminosity give consistent results, as can be seen 
for Ha in Figure 23 (an almost identical comparison results for [On]). This indicates the high level of self-consistency 
between the photometric and spectroscopic data in the SDSS, and suggests that both the measurements and the methods 
are self-consistent. There remains a small systematic discrepancy, though, such that SFRs calculated using the aperture 
corrections based on Equation Al are about 15% larger. This is the case for both Ha and [On], and may be related to 
the spectrophotomctric calibration. The SDSS takes spectra during conditions which are not deemed "photometric" . The 
seeing conditions are thus typically worse for the spectroscopic data than for the photometric data. Early versions of the 
photometric pipeline software, used to photometrically calibrate the spectra, did not take these seeing differences into 
account. In the DR1, however, all photometric data used in calibrating the spectra are convolved to 2" seeing (typical 
of the seeing conditions for spectroscopic observations). The currently available spectra have, unfortunately, not yet 
been recalibrated, resulting in spectrophotometric magnitudes fractionally brighter than what one would expect from the 
photometrically measured fiber magnitudes (by « 0.1 magnitude in the mean). The SFRs using the aperture correction of 
Equation A2 and the emission line fluxes (Equation B3) thus slightly overestimates the correction (since the line fluxes are 
slightly overestimated). The alternative SFR estimate (Equation B2), which uses the line EWs and absolute magnitudes, 
is insensitive to this issue. 

We investigated the aperture corrections of Equation Al using Petrosian and fiber magnitudes in both r and z. The 
r magnitudes were seen to give a distribution with less scatter than the z magnitudes, and as a result we chose to apply 
those in preference. For the estimation of [On] luminosities, the u-band Petrosian and fiber magnitudes were used in 
applying Equation Al, but very little difference results if the r-band magnitudes are retained instead. 

To emphasise the extent of the aperture corrections we show how they vary with several parameters. The logarithm of 
the aperture corrections is shown as a function of galaxy size (the Petrosian radius) in Figure 24, and Figure 25 shows 
the aperture correction as a function of redshift and of SFRi. 4 gh z - It can be seen from these Figures that the aperture 
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corrections are typically at least a factor of two and can be as much as an order of magnitude. Finally, the variation in 
the ratio of the SFRs from Ha and 1.4 GHz luminosities with the aperture correction is shown in Figure 26. While the 
relation is almost flat, as would be desired for an aperture correction introducing no bias in the derived SFR, there is a 
measurable positive slope, which reflects the implicit assumption of a uniform SF distribution throughout each galaxy. In 
systems with the largest aperture corrections, the SFRna is slightly overestimated, as already seen in Figure 10. There 
are a small number of systems with very small aperture corrections, and another small population at all values of aperture 
correction, where SFRi. 4 gh z appears to be strongly overestimated, possibly due to the presence of low luminosity AGN 
components contributing to the radio emission. Apart from these systems, at low aperture corrections the median ratio 
of the two SFRs is unity. 

B. SDSS SFR FORMULAE 

For ease of reference, the formulae for deriving SFRs using the measured SDSS parameters, and the SFR calibrations 
used, are all collected together here. The sections in which each formula is derived are also given. All SFRs are calibrated 
based on a Salpeter IMF and a mass range of 0.1 < M© < 100. 
The Ha luminosity to SFR calibration used is 

SFR Ha (M 0yr - 1)= i27 ^ O 3 4w . (Bl) 

For Ha SFRs, using the aperture correction method of Equation 5, the derivation in §3.2.2 gives 

SFR Hq (M© yr" 1 ) = (EW(Ha) + EW C ) Kr°- 4 ( M '- 34 - 10 ) 

3 x 10 18 ( Sua/SuA 2114 1 (B2) 
(6564.61(1 + z)) 2 \ 2.86 / 1.27 xlO 34 ' 1 ' 

where Sn a and Sup are the stellar absorption corrected line fluxes, calculated as in Equation 4. The exponent on the 
Balmer decrement term (in all the equations given here) is equal to fc(A)/[fc(H/3) — fc(Ha)], and depends on the assumed 
obscuration curve. For obscuration corrections to emission line luminosities, we assume the obscuration curve of Cardelli 
et al. (1989) as recommended by Calzetti (2001). EW C = 1.3 A is a reasonable approximation for the stellar absorption 
correction when using the SDSS pipeline spectral line measurements, and corresponds roughly to a 2.6 A EW stellar 
absorption in the SF galaxies. Using the alternative aperture correction given in Appendix A results in 

SFR HQ (M© yr" 1 ) = 4nD? S Ha 10^ 4 (— -> (^f^) ^ L27 ' 1Q34 , (B3) 

where Di is the luminosity distance, and Su a is the stellar absorption corrected Ha line flux. 
The [On] luminosity to SFR calibration used is 

SFRroiii (M Q yr- 1 ) = L[OIJl „ — , (B4) 

lonj v oj ) 2.97 xlO 33 W v ' 

where £[oii] incorporates the obscuration correction valid for the wavelength of Ha. For [On] SFRs the derivation of § 3.3 

gives 

SFR [OII] (M Q yr- 1 ) = EW(OII) i(r - 4 ( M "- 34 - 10 ) 

3 xlO 18 {SnJSnpY 114 1 (B5) 



(3728.30(1 + z)) 2 \ 2.86 J 2.97 xlO 33 ' 
and using the alternative aperture correction given in Appendix A results in 

SFR [OII] (M© yr -1 ) = A,D 2 F [OII] lO" ^"") (^f^)" 14 2 . 97 ' 10 3 3 - ( B6 ) 
The w-band luminosity to SFR calibration used is 



1.186 



SFRu (Moyr- 1 ) = ( ^ r ) . (B7) 

U V oy ' 1 1.81 x 10 21 WHz- 1 ' 



The derivation given in §3.5 gives 



2.061 \ 1186 



SFR (M vr-M - f 10-°- 4(M "- 3410) fW, , 
SFRu (M© yr ) - ^ ^ —— j j . (B8) 

The exponent on the Balmer decrement term here uses the obscuration curve of Calzetti et al. (2000), and incorporates 
the factor of 0.44 necessary for obscuration corrections of the stellar continuum (see also Calzetti 2001). 
For completeness, the SFR calibrations from 1.4 GHz and FIR luminosities that give consistent SFR estimates with the 
above formulae are also given here (from Bell 2003). The calibration for 1.4 GHz luminosities is 

( M vr-ls- j £l.4GHz/[l-81 X lO^tWHz" 1 )] L 1AGRz > L c . . 

SFRi.4GHz {Mq yr )-| Xl 4GHz /[( .l + OJ^l^/Lc) 03 ) 1.81 x 10 21 (W Hz" 1 )] L 1AGHz < L c , (B9) 
with L c = 6.4 x 10 21 WHz -1 , and that for FIR luminosities is 

SFR Fm (M q yr -1 ) - L Fm (l + ^2.186 x 10 35 (W)/L F ir)/[1.85 x 10 36 (W)]. (B10) 
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Table 1 

Scatter of SFR indicators relative to SFRi. 4 gh z 



Indicator 


rms scatter 


Notes 


FIR 


±40% 


Using calibration from Equation 7 


Ha 


±60% 


Using obscuration curve from Cardelli et al. (1989) 


[Oil] 


±70% 


Using F OI1 /F Ha = 0.23 


u-band 


±70% 


Using derived calibration of Equation B7 
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0: SF 
1: AGN 

2: Unclassified 


2 2 2 
638 


Emission: 2407 
Quiescent: 672 
Total: 3079 
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Fig. 1. — Distribution of the spectroscopic diagnostic classifications for all 3079 DR1 main galaxies with 1.4 GHz FIRST detections. The 
three numbers at the top of each box are the flags (0, 1, or 2) indicating the classification in the three spectral diagnostic diagrams. Below 
these is the number of sources with each particular combination of flags (for example, there are 791 sources with for each flag). There are 
672 "quiescent" systems, with none of the necessary emission lines for classification. 
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Fig. 2. — Redshift distribution of the 1.4 GHz FIRST detections in the sample. The histograms indicate the proportions classified as SF, 
AGN and "quiescent." This is done in a cumulative fashion, adding the histogram for each population onto the previous total, to emphasise 
the relative proportions in each bin. 
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Fig. 3. — 1.4 GHz flux density distribution for the FIRST detections in the sample. As in the previous Figure, the histograms arc cumulative, 
and indicate the proportions classified as SF, AGN and "quiescent." 
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Fig. 5. — Galaxy size given by the Petrosian radius in r-band as a function of SFRi.4ghz- Here it can be explicitly seen that the FIRST 
SFRs below about 10 Mq yr _1 belong to progressively larger galaxies, implying that these SFRs are progressively underestimated. 
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Fig. 6. — (a) Ratio of NVSS to FIRST luminosities for SF galaxies as a function of Petrosian radius in r-band. This shows that the 
majority of sources smaller than about r = 10" have comparable luminosities (and flux densities) from NVSS and FIRST, although NVSS 
still appears to be fractionally larger (by about 5% to 10%) on average. Above 10", though, the NVSS flux density is consistently higher than 
that from FIRST. There are a handful of sources with r < 10" with significantly higher NVSS flux densities, and these are likely to be cases 
where one or more nearby galaxies have entered the NVSS beam in addition to the galaxy for which the FIRST flux density is measured, (b) 
Comparison of the SFRs derived using NVSS and FIRST 1.4 GHz luminosities, for all SF galaxies having r > 10". The greater sensitivity 
of the NVSS to extended radio emission is clearly seen in the comparison. The error bars in the upper left indicate the typical uncertainty 
in the measurements. They represent random error only, no systematic errors (such as the uncertainty in the SFR calibration) are included. 
The lower panel shows the residuals from the one-to-one line. 
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Stellar absorption EW A 

Fig. 7. — Histogram of intrinsic stellar absorption EWs at Ha for the complete sample of SF galaxies. The median value is 2.6 A. Inset: 
Sample SDSS spectrum indicating the effect of the pipeline measurement of the H/3 emission line, illustrating the incomplete extent to which 
stellar absorption reduces the measured emission. The points are the actual measured SDSS spectrum. The solid line is the SDSS pipeline 
Gaussian fit to the emission, and the dashed lines are a double Gaussian fit, made by fixing an absorption component to have an EW= 3.0 A 
the same EW as the absorption measured at H5 in this spectrum. The black shaded region indicates the true flux missed by the pipeline 
Gaussian fit in the measurement of the emission, while the grey shaded regions indicate excess flux added by the pipeline Gaussian fit. The 
black area minus the grey area is the true amount by which the pipeline fit underestimates the emission. It can clearly be seen that this area 
is smaller (by at least a factor of two) than the area of the stellar absorption component. Given a median stellar absorption of EW= 2.6 A 
at Ho, the SDSS pipeline measurements, typically, are diminished by only about EW= 1.3 A. 
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Fig. 8. — Balmer decrements, corrected for stellar absorption as described in the text, as a function of SFRi.4ghz- The different symbols 
indicate ranges in Ha equivalent width, as shown. The error bars in the upper left indicate the typical uncertainty in the measurements. They 
represent random error only, no systematic errors (such as the uncertainty in the SFR calibration) are included. The dashed line indicates 
the 2.86 value expected from case B recombination (Brocklchurst 1971), the solid line indicates the Balmer decrements predicted from the 
SFR-dcpcndcnt obscuration derived by Afonso ct al. (2003), and the dot-dashed line those from Hopkins ct al. (2001). (These relations have 
been converted to be consistent with the 1.4 GHz SFR calibration and cosmology used here.) While the empirical correction of Hopkins et 
al. (2001) is clearly an underestimate for the sample on average (although it may be reasonable for systems with EW(Hq) > 70 A), that of 
Afonso ct al. (2003) seems to be somewhat of an overestimate, on average, for this sample. This is likely to reflect the radio-selected nature 
of the sample from which this relation was derived, suggesting that the present sample may be missing a number of more highly obscured 
systems. 
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Fig. 9. — SFRs from Ha luminosity as a function of SFRs from FIRST 1.4 GHz luminosity. Only galaxies spectroscopically classified as SF 
are shown. All panels show the 771 SF galaxies remaining when NVSS measurements are used for galaxies having r > 10". Panel (a) shows 
SFRh q uncorrected for aperture or obscuration effects. Panel (b) shows the aperture corrected SFRh q prior to the obscuration correction, 
calculated from Equation 5, and panel (c) adds the obscuration correction using the Balmer decrement to give the fully corrected SFRhq,. 
Panel (d) uses the method of Afonso ct al. (2003) rather than the Balmer decrement for making the obscuration correction, for comparison. 
The rms deviation either side of the one-to-one line in (c) is 0.21 dex, a factor of 1.6. The error bars in the upper left of panel (c) indicate 
the typical uncertainty in the measurements. Again, they include random error only. 
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Fig. 10. — Ratio of SFRs from 1.4 GHz and Ha luminosities as a function of redshift. The horizontal line shows the one-to-one relationship, 
and the curved line indicates the effect of the flux density limit of the FIRST survey detections (about 0.75 mjy) combined with the approximate 
upper limit of the measured Ha SFRs (about 100 Mq yr _1 ). The symbol sizes reflect the apparent size of the object, based on the Petrosian 
radius as indicated in the Figure. The error bar in the lower right indicates the typical uncertainty in the SFR ratio. This is again random 
error from the measurements only. 
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Fig. 11. — SFR[on] compared with SFRi.4GHz and SFRjja- The rms deviation either side of the one-to-one line is 0.22 dex, a factor of 
1.7, for the comparison with SFRi.4gjjz> an d is 0.15 dex, a factor of 1.4, for the comparison with SFRh q . Interestingly, the systems with 
SFR[on] < lMQyr -1 seem to all show significant overestimates in SFRi.4qhz (while the SFRjjq is consistent with the SFR[on])- These 
systems are likely to be hosting a heavily obscured AGN, which dominates the radio emission but is not detectable through the optical 
spectroscopic signature. The error bars in the upper left of both panels indicate the typical random error in the measurements. 
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Fig. 12. — SFRfir from Equation 7 compared with (a) SFRi.4GHz and (b) SFRhq. SFRjja is calculated as given in Equation B2. The 
rms deviation either side of the one-to-one line is 0.15 dex, a factor of 1.4, for the comparison with SFRi.4ghz, and 0.22 dcx, a factor of 1.7 
for the comparison with SFRh q ■ The error bars in the upper left of both panels indicate the typical random error in the measurements. It is 
interesting to note in panel (b) that for SFRpjR < 7M0yr _1 , almost all the points show relatively high SFRjia/SFRpiR. These galaxies are 
mostly nearby systems with large aperture corrections for the Ha SFR estimate, and it is likely that the SFRjja is somewhat overestimated 
for these objects. 
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Fig. 13. — (a) Comparison between the M-band luminosities and the Ho SFRs, for the complete sample of 2625 SF galaxies in the DR1. 
Luminosities have been corrected for obscuration using the stellar absorption corrected Balmer decrements. The SFRjj a estimates have also 
been corrected for aperture effects. The solid line shows the ordinary least squares bisector fit described in the text, (b) The resulting SFR 
comparison after applying the new SFRtj calibration. The rms deviation either side of the one-to-one line is 0.13 dex, a factor of 1.3. The 
error bars in the upper left of both panels indicate the typical random error in the measurements. 
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Fig. 14. — Comparison between the «-band and Ha SFRs, for all 21649 SF galaxies with the necessary measurements in the DR1. Despite 
the incompleteness when the entire sample is considered, the derived calibration is very similar to the ordinary least squares bisector line fit. 
There can still be seen, nevertheless, as a slight overall apparent bias to overestimated SFRtj at high SFRs, although the two SFR estimates 
are in fact consistent (see Figure 13). The error bars in the upper left indicate the typical random error in the measurements. 
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Fig. 15. — SFRs from u-band luminosity compared with the SFRs from FIRST 1.4 GHz luminosity. Shown are the 771 SF galaxies remaining 
after NVSS measurements are used to replace those from FIRST for galaxies with r > 10". The upper panel (a) shows SFRu before the 
obscuration correction is applied. The lower panels incorporate the obscuration correction to the SFRu, directly from the (stellar absorption 
corrected) Balmer decrement in (b), and using the method of Afonso et al. (2003) in (c). The error bars in the upper left of panel (b) indicate 
the typical random error in the measurements. The slight apparent offset in (b) towards lower SFRtj estimates is a result partially of the 
incompleteness of this sample although the same effect is present, to a lesser extent, if the radio detected galaxies arc restricted to those in 
the complete sample (see discussion in §4.2, and Figure 18). 
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Fig. 16. — Histogram showing the range of absolute magnitude, M g , spanned for the complete sample of 2625 SF galaxies. The radio 
detected galaxies within the complete sample are indicated by the filled histogram (note that different bin sizes have used, to emphasise the 
shape of the distribution for the radio detections). For the whole complete sample the median M g is —20.6, while for the radio detected 
galaxies it is —20.8, although the distribution for the radio detected systems is much broader. Gaussian fits to the two distributions here 
gives a FWHM of 0.9 mag for the complete sample, and 1.3 mag for the radio detected galaxies. 
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Fig. 17. — Histograms showing, for the complete sample of 2625 SF galaxies, the difference in (a, b) colors, (c) inverse of the concentration 
index, 1/c, and (d) D4ooo- The radio detected galaxies within the complete sample are indicated by the filled histograms. The dashed 
histograms in each panel show the distribution for the complete sample of all galaxies, 24444 objects, not restricted to star-forming systems. 
The heights of these histograms have been scaled down by a factor of ten to emphasise the relative shapes of the distributions, the bimodal 
nature of which is clear, being split clearly into red and blue populations, with the blue population hosting the star-forming systems. It can 
be seen that the radio detected SF population has redder colors, smaller 1/c and larger D4000 values than the SF population as a whole. This 
suggests are larger contribution to the optical emission, on average, from old-stellar populations in radio detected SF systems. 
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Fig. 18. — Comparison between SFRtj and SFRjja for the complete sample (dots), indicating the location of the radio detected systems 
(filled circles). The error bars in the upper left indicate the typical random error in the measurements. It can be seen that, while high SFR 
systems (ss lOOM0yr — 1 ) seem to be uniformly detected at 1.4 GHz, for moderate SFR galaxies (ss lOMoyr -1 ) radio detection seems to 
preferentially select the lower-luminosity n-band systems. 
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Fig. 19. — Histogram showing the (emission line) obscuration in magnitudes at the wavelength of Ha for the complete sample of SF 
galaxies. The radio detected galaxies are shown as the filled histogram. The median obscuration for the complete optically selected sample is 
1.2 magnitudes, while for the radio detected galaxies, it is 1.6 magnitudes. Gaussian fits to these histograms give FWHM values of 1.0 mag 
for the complete sample, and 1.2 mag for the radio detected galaxies. 
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Fig. 20. — (a) Histogram of F[ou]/Fna for the full sample of 22577 SF galaxies from DR1. The 752 radio detected galaxies in the full 
sample are shown as the filled histogram. The median value for the whole (incomplete) sample is -F[on]/^Ha = 0.38, while for the radio 
detected sources it is 0.23. (b) The distribution for the complete sample of SF galaxies. The median value for the complete, optically selected 
sample is 0.26. The radio detected galaxies in the complete sample are shown as the filled histogram, with a median of 0.21. 
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Fig. 21. — Ratio of [On] to Ha line fluxes as a function of EW(Ho) for the complete sample of SF galaxies. A weak trend can be seen 
for higher flux ratios to be present in systems of higher EW. The error bars in the upper left indicate the typical random error in the 
measurements. 
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Fig. 22. — Relationship between the F[qh] /-Fka ratio and the Balmer decrement for the complete sample of SF galaxies. The trend seen 
is to be expected, given that the -Ffoiq/^Ha flux ratio is constructed prior to obscuration corrections, and hence will trace the obscuration 
depending on the obscuration curve used. The error bars in the upper right indicate the typical random error in the measurements. The 
three lines shown in (a) indicate the trend expected from the obscuration curve of Cardelli et al. (1989) for intrinsic S^qu]/ Sn a flux ratios of 
0.3, 1.0 and 2.0, from bottom to top respectively. 
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Fig. 23. — Comparison of the H« SFRs derived using the different methods described in Appendix A. The abscissa shows SFR calculated 
using Equation B2, based on EW(Hq) and absolute r-band magnitude, and the ordinate shows SFR calculated using Equation B3, based on 
apparent r-band magnitude and fiber magnitude. Both SFR estimates have been corrected for obscuration using the Balmer decrement. The 
error bars in the upper left indicate the typical random error in the measurements. The lower panel shows the residuals from the one-to-one 
line. 




Fig. 24. — The aperture correction (from Equation Al) applied to the fiber-based SFR estimates shown as a function of galaxy size as given 
by the Petrosian radius. Clearly the larger galaxies require a larger aperture correction. 
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Fig. 25. — The aperture correction (from Equation Al) as a function of (a) redshift and (b) SFRi.4GHz- The largest galaxies with greatest 
aperture corrections lie at the lowest redshifts, and have the lowest radio derived SFRs. It can be inferred that below about 10 Mq yr~ x the 
FIRST SFRs become progressively further underestimated. This is shown explicitly in Figure 5. 
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Fig. 26. — The ratio of SFRs from Ha and 1.4 GHz luminosities as a function of the aperture correction (from Equation Al). The implicit 
assumption of a uniform SF distribution made through the aperture correction results in the slight positive slope seen in this relation. 



